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Taking the examples of Legendre and Hermite orthogonal polynomials, we show 
how to interpret the fact that these orthogonal polynomials are moments of other 
orthogonal polynomials in terms of their associated Riordan arrays. We use these 
means to calculate the Hankel transforms of the associated polynomial sequences. 

1 Introduction 



In this note, we shall re-interpret some of the results of Ismail and Stanton [15, 16] in terms 
of Riordan arrays. These authors give functional [15] whose moments are the Hermite, La- 
' guerre, and various Meixner families of polynomials. In this note, we shall confine ourselves 

to Legendre and Hermite polynomials. Indeed, the types of orthogonal polynomials repre- 
sentable with Riordan arrays is very limited (see below), but it is nevertheless instructive 



Q ! to show that a number of them can be exhibited as moments, again using (parameterized) 

Riordan arrays. 

The essence of the paper is to show that a Riordan array L (either ordinary or exponential) 
defines a family of orthogonal polynomials (via its inverse L~^) if and only if its production 
><1 : matrix [8, 9, 10] is tri-diagonal. The sequence of moments /i„ associated to the family of 

, orthogonal polynomials then appears as the elements of the first column of L. In terms of 

generating functions, this means that if L = {g, f) (or L = [g, /]), then g{x) is the generating 
function of the moment sequence. By defining suitable parameterized Riordan arrays, we 
can exhibit the Legendre and Hermite polynomials as such moment sequences. 

While partly expository in nature, the note assumes a certain familiarity with integer 
sequences, generating functions, orthogonal polynomials [5, 12, 30], Riordan arrays [25, 29], 
production matrices [10, 23], and the Hankel transform of sequences [2, 7, 20]. We provide 
background material in this note to give a hopefully coherent narrative. Many interesting 
examples of sequences and Riordan arrays can be found in Neil Sloane's On-Line Encyclo- 
pedia of Integer Sequences (OEIS), [27, 28]. Sequences are frequently referred to by their 
OEIS number. For instance, the binomial matrix B ("Pascal's triangle") is A007318 . 

The plan of the paper is as follows: 
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2 Preliminaries on integer sequences and Riordan ar- 
rays 

For an integer sequence a„, that is, an element of Z^, the power series /(x) = ^^q*^"-^" 
is called the ordinary generating function or g.f. of the sequence, is thus the coefficient 
of x" in this series. We denote this by a„ = For instance, Fn = [x^]-[z§z^ is 

the n-th Fibonacci number A000045 , while C„ = [x"] ^"^^ -^ is the n-th Catalan number 
A000108 . The article [21] gives examples of the use of the operator [x"]. We use the notation 

^ [3,n]^ ^Yie sequence 1, 0, 0, 0, ... , A000007 . Thus 0" = [n = 0] = 5^,0 = 0- Here, 
we have used the Iverson bracket notation [13], defined by \P] = 1 if the proposition V is 
true, and \P] = if P is false. 

For a power series /(x) = ^^o*^"^" with /(O) = we define the reversion or composi- 
tional inverse of / to be the power series /(x) such that /(/(x)) = x. We shall sometimes 
write this as / = Rev/. 

For a lower triangular matrix {an,k)n,k>Q the row sums give the sequence with general term 
^^=0 ^n^k while the diagonal sums form the sequence with general term 



The Riordan group [25, 29], is a set of infinite lower-triangular integer matrices, where 
each matrix is defined by a pair of generating functions g{x) = 1 + gix + g2x'^ + ■ ■ ■ and 
/(x) = /ix + /2X^ + ■ ■ ■ where /i 7^ [29]. We assume in addition that /i = 1 in what 




fc=0 
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follows. The associated matrix is the matrix whose i-th column is generated by g{x)f{xY 
(the first column being indexed by 0). The matrix corresponding to the pair f is denoted 
by {g, f) or 7l{g, f). The group law is then given by 

{gJ)-{hJ) = {gJ){hJ) = {g{hof)Jof). 

The identity for this law is / = {l,x) and the inverse of {g,f) is {g,f)~^ = (l/(fi' o /),/) 
where / is the compositional inverse of /. 

A Riordan array of the form {g{x),x), where g{x) is the generating function of the se- 
quence ttn, is called the sequence array of the sequence a„. Its general term is a„_fc (or more 
precisely, [k < n]a„_fc). Such arrays are also called Appell arrays as they form the elements 
of the so-called Appell subgroup. 

If M is the matrix ((?,/), and a = (oq, ai, . . .)' is an integer sequence with ordinary gener- 
ating function A (x), then the sequence Ma has ordinary generating function g{x)A{f{x)). 
The (infinite) matrix {g, f) can thus be considered to act on the ring of integer sequences 
by multiplication, where a sequence is regarded as a (infinite) column vector. We can 
extend this action to the ring of power series Zi[[x]] by 

{g, /) : A{x) ^ {g, /) ■ A{x) = g{x)A{f{x)). 

Example 1. The so-called binomial matrix B is the element {jz^, yz^) of the Riordan 
group. It has general element (^), and hence as an array coincides with Pascal's triangle. 
More generally, B"* is the element ( iJ^^ , j^^) of the Riordan group, with general term 
(^)m"~*^. It is easy to show that the inverse B""^ of B"* is given by ( i^^nx ' i+mx ) • 

Example 2. If a„ has generating function g{x), then the generating function of the sequence 

LfJ 

bn = an-2k 
k=0 

is equal to 

while the generating function of the sequence 



is equal to 
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The row sums of the matrix {g, f) have generating function 



(9 J) 



X 



9[x) 
1 - fix) 



while the diagonal sums of {g, /) (sums of left-to-right diagonals in the North East direc- 
tion) have generating function g{x)/{l — xf{x)). These coincide with the row sums of the 
"generalized" Riordan array {g,xf): 



{9,xf) 



9[x) 



X 



xf{x)' 



For instance the Fibonacci numbers -F„+i are the diagonal sums of the binomial matrix B 
given by 
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while they are the row sums of the "generalized" or "stretched" [6] Riordan array ( , -^-^ j : 



/ 1 ... \ 

1 

1 1 

1 2 

13 10 

1 4 3 



V 



Each Riordan array {g{x), f{x)) has bi-variate generating function given by 

9{x) 

1 - yf{x)' 

For instance, the binomial matrix B has generating function 

1 



1 

l-x 



l-yj^ l-x{l+yy 



For a sequence ao,ai,a2, . . . with g.f. g{x), the "aeration" of the sequence is the sequence 
ao, 0, ai, 0, a2, . . . with interpolated zeros. Its g.f. is gi^x"^)- 
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The aeration of a (lower-triangular) matrix M with general term rriij is the matrix whose 
general term is given by 

1 + (-1)'-^' 



z+j i — j 



where , is the i, j-th element of the reversal of M: 



m 



In the case of a Riordan array (or indeed any lower triangular array), the row sums of the 
aeration are equal to the diagonal sums of the reversal of the original matrix. 

Example 3. The Riordan array (c(x^), xc(x^)) is the aeration of {c{x) , xc{x)) A033184 . Here 

1 - VI -4x 



c{x) 



2x 



is the g.f. of the Catalan numbers. Indeed, the reversal of {c{x),xc{x)) is the matrix with 
general element 

+ k\n — k + 1 
k ) n+1 ' 



[k<n+l] 



which begins 



/ 1 

1 1 

1 2 2 

1 3 5 5 

1 4 9 14 14 

1 5 14 28 42 42 



V : : : : : : 

This is A009766 . Then (c(x^), xc(x^)) has general element 



1\ A; + ll + (-l) 



n 



^ Jn + 1 



n—k 



and begins 



/ 1 

1 

10 10 

2 1 

2 3 1 

5 4 1 
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This is A053121. Note that 



(c(x^), a;c(x^)) 



X 



1 + 1 + 
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We observe that the diagonal sums of the reverse of (c(x), xc(x)) coincide with the row sums 
of (c(x^), xc(x^)), and are equal to the central binomial coefficients (^j) A001405 . 

An important feature of Riordan arrays is that they have a number of sequence charac- 
terizations [4, 14]. The simplest of these is as follows. 

Proposition 4. [14, Theorem 2.1, Theorem 2.2] Let D = [dn^k] be an infinite triangular 
matrix. Then D is a Riordan array if and only if there exist two sequences A = [a^, Oi, a2, . . .] 
and Z = [zq, Zi, Z2, ■ ■ ■] with 7^ 0, 7^ such that 

• dn+l,k+l = 0'jdn,k+j, (k, n = 0, 1, . . .) 

• dn+lfi = YlJLo ^j'^n,j^ (n = 0, 1, . . .). 

The coefficients ao,ai,a2, . . . and zo,zi,Z2,... are called the A-sequence and the Z- 
sequence of the Riordan array D = {g{x), f{x)), respectively. Letting A{x) be the generating 
function of the A-sequence and Z{x) be the generating function of the Z-sequence, we have 



3 Orthogonal polynomials and Riordan arrays 

By an orthogonal polynomial sequence {pn{x))n>o we shall understand [5, 12] an infinite 
sequence of polynomials Pnix), n > 0, of degree n, with real coefficients (often integer 
coefficients) that are mutually orthogonal on an interval [xq, xi] (where xq = — c>o is allowed, 
as well as xi = oo), with respect to a weight function w : [xq, xi] — )■ M : 



Xo 



where 



Pn{x)Pm{x)w{x)dx = 5nm\/ KK 



x\ 

p'^{x)w{x)dx = hn- 



' xo 

We assume that w is strictly positive on the interval (xo,xi). Every such sequence obeys a 
so-called "three-term recurrence" : 

Pn+l{x) = (a„X + bn)Pn{x) - C„p„_i(x) 

for coefficients a„, 6„ and c„ that depend on n but not x. We note that if 

Pj{x) = kjx'-' + k'jX^~^ + . . . j = 0, 1, . . . 

then 



Since the degree of Pn{x) is n, the coefficient array of the polynomials is a lower triangular 
(infinite) matrix. In the case of monic orthogonal polynomials (where kn = 1 for all n) the 
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diagonal elements of this array will all be 1. In this case, we can write the three-term 
recurrence as 

Pn+lix) = {X- a„)Pnix) - /3„p„_l(x), Po{x) = 1, = X - Oq. 

The moments associated to the orthogonal polynomial sequence are the numbers 

Jxo 

We can find Pn{x), an and /3„ (and in the right circumstances, w{x) - see the Appendix) from 
a knowledge of these moments. To do this, we let A„ be the Hankel determinant |/^j+j|"j>o 
and An^x be the same determinant, but with the last row equal to 1, x, x^, . . .. Then 



term Hu,+v, ■ Let 



Vl ... Vk 



More generally, we let ( ' ' ' ) be the determinant of Hankel type with («, j)-th 



. ^^/01...n\ / 01...n — 1 n 

01. ..n/ V01...n — ln + 1 



Then we have 



" A„ A„.i' A2 • 



We shall say that a family of polynomials {Pn(2;)}n>o is [31] formally orthogonal, if there 
exists a linear functional £ on polynomials such that 

1. Pnix) is a polynomial of degree n, 

2. C{pn{x)Prn{x)) = for m 7^ n, 

3. ^ 0. 

Consequences of this definition include [31] that 

£(x™p„(x)) = Kn^mn, < m < U, ft„ 7^ 0, 

and if q{x) = Y.l=o"'kPk{x), then Uk = C{qpk)/C{pl). 

The sequence of numbers /i„ = C{x"') is called the sequence of moments of the family 
of orthogonal polynomials defined by C Note that where a suitable weight function w{x) 
exists, then we can realize the functional C as 

C{p{x)) = / p{x)w{x) dx. 
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If the family Pn{x) is an orthogonal family for the functional C, then it is also orthogonal 
for cC, where c 7^ 0. In the sequel, we shall always assume that /io = C{po{x)) = 1. 

The following well-known results (the first is the well-known "Favard's Theorem"), which 
we essentially reproduce from [18], specify the links between orthogonal polynomials, three 
term recurrences, and the recurrence coefficients and the g.f. of the moment sequence of the 
orthogonal polynomials. 

Theorem 5. [18] (Cf. [31], Theoreme 9 on p. 1-4, or [32], Theorem 50. i;. Let (p„(x))„>o 
he a sequence of monic polynomials, the polynomial Pn{x) having degree n = 0, 1, . . . Then 
the sequence is (formally) orthogonal if and only if there exist sequences (a„)„>o and 

(/3n)n>i with /3n 7^ for all n > 1, such that the three-term recurrence 

Pn+l = (X - an)Pn{x) " /3nPn-l(x), for n > 1, 

holds, with initial conditions po{x) = 1 and pi{x) = x — a^. 

Theorem 6. [18] (Cf [31], Proposition 1, (7), on p. V-5, or [32], Theorem 51. i;. Let 
{Pn{x))n>Q be a sequence of monic polynomials, which is orthogonal with respect to some 
functional C. Let 

Pn+i = {x- an)Pn{x) - /3„p„-i(x), for > 1, 

he the corresponding three-term recurrence which is guaranteed hy Favard's theorem. Then 
the generating function 

oo 

9{x) = ^^ikX^ 

for the moments fik = C{x^) satisfies 

9{x) = 



1 — aox 



P2X 

1 — aix 



[^3X' 

1 — a2X 



1 — — ■ ■ ■ 
Given a family of monic orthogonal polynomials 

Pn+l{x) = (X- an)Pn{x) " (3nPn-l{x), Po{x) = 1, Pi{x) = X - Oq, 

we can write 



Then we have 



Pnix) = '^an,kx''. 

k=0 

n+1 n n—1 

an+i,kx'' = {x - an) ^ an,kx'' - f3n^ a,n-l,kX^ 

k=0 k=0 k=0 
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from which we deduce 

C^n+l.O = —OinO'nfi ~ /^n^n-l.O (2) 

and 

an+l,fc = fln.fc-l — OtnCtn,k — I^n0'n-l,k (3) 

We note that if a„ and /3„ are constant, equal to a and /3, respectively, then the sequence 
(1, —a, — /3, 0,0,.. .) forms an A-sequence for the coefficient array. The question immediately 
arises as to the conditions under which a Riordan array [g, f) can be the coefficient array of 
a family of orthogonal polynomials. A partial answer is given by the following proposition. 

Proposition 7. Every Riordan array of the form 

1 X 



1 + rx + sx"^ 1 + rx + sx"^ ^ 
is the coefficient array of a family of monic orthogonal polynomials. 
Proof. The array ( i+ , i+^x+sx^ ) t^'^] ^ C-sequence C{x) = Y.n>o'^nX'^ given by 

X X 

1 + rx + sx"^ 1 — xC{x) ' 

and thus 

C{x) = — r — sx. 

This means that the Riordan array f-n — \ — 2, t, — — 2) is determined by the fact that 

\ l+rx+sx^ ' 1+rx+sx'^ I 

a„+i,fc = a„,fc-i + ^ CiO'n-i^k for n, /c = 0, 1, 2, . . . 

i>0 

where a„ _i = 0. In the case of f-n — \ — 2 , -r^ — -^ — 2) we have 

On+l,fc = CLn,k-l ~ ^0'n,k ~ San-l,k- 

Working backwards, this now ensures that 

Pn+l{x) = (X - r)pn{x) - SPn^l{x), 

where = X]fc=o '^"■fc-^"- "^^^ result now follows from Theorem 5. □ 

We note that in this case the three-term recurrence coefficients a„ and /3,„ are constants. 
We have in fact the following proposition (see the next section for information on the Cheby- 
shev polynomials). 

Proposition 8. The Riordan array (^q:;:;^:]:^, Y+T^jii^) ^■^ coefficient array of the mod- 
ified Chebyshev polynomials of the second kind given by 



Pn{x) = {V~srUn (^) ' n = 0,1,2, 



9 



Proof. The production array of 



1 



) ^ is given by 



/ r 1 

s r 1 

s r 1 

s r 1 

s r 1 

s r 



V : : : : : : •• / 

The result is now a consequence of the article [11] by Elouafi, for instance. 



□ 



The complete answer can be found by considering the associated production matrix of a 
Riordan arrray, in the following sense. 

The concept of a production matrix [8, 9, 10] is a general one, but for this work we find 
it convenient to review it in the context of Riordan arrays. Thus let P be an infinite matrix 
(most often it will have integer entries). Letting fq be the row vector 



:i, 0,0,0,...) 



we define rj = rj_iP, i > 1. Stacking these rows leads to another infinite matrix which we 
denote by Ap. Then P is said to be the production matrix for Ap. 
If we let 



then we have 



M^ = (1,0,0,0,...,0,...) 



u^P 

Ap = ^Tp2 



\ 



J 



and 



lAp = ApP 

where I = {Si+i,j)i,j>o (where S is the usual Kronecker symbol): 

/ 1 ... \ 

1 

1 

1 

1 





V 



We have 



P = Ap^IAp. 



(4) 
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Writing Ap = lAp, we can write this equation as 



P = Ap'Ap. 



(5) 



Note that is a "beheaded" version of Ap; that is, it is Ap with the first row removed. 

The production matrix P is sometimes [23, 26] called the Stieltjes matrix Sap associated 
to Ap. Other examples of the use of production matrices can be found in [1], for instance. 
The sequence formed by the row sums of often has combinatorial significance and is 
called the sequence associated to P. Its general term a„ is given by = u^P"'e where 

/ 1 \ 
1 
1 

V '■■ I 

In the context of Riordan arrays, the production matrix associated to a proper Riordan array 
takes on a special form : 

Proposition 9. [10, Proposition 3.1] Let P he an infinite production matrix and let Ap he 
the matrix induced hy P. Then Ap is an (ordinary) Riordan matrix if and only if P is of 
the form 

\ 





«0 
















«! 


«o 













«2 


«! 


«o 










"3 


^2 


a\ 


«o 





^4 


04 


"3 


a2 


«! 






"5 


04 


"3 


«2 





V 



where 7^ 0, «o 7^ 0. Moreover, columns and 1 of the matrix P are the Z- and A- 
sequences, respectively, of the Riordan array Ap. 

We recall that we have 



A{x) 



X 



Zix) 



fix)' fix 
Example 10. We consider the Riordan array L where 



gifix)) 



Xx 



jJX 



X 



1 + ax + hx"^ ' 1 + ax + hx"^ 
The production matrix (Stieltjes matrix) of 



1 - Xx 



fix 



X 



1 + ax + hx"^ 1 + ax + hx'^ 



11 



is given by 



P = Sr 



( a+A 10 

6 + a 1 

6 a 1 

6 a 1 

6 a 1 










6 



a 



We note that since 



Xx 



fix 



X 



1 + ax + bx'^ ' 1 
(1 — Ax — jix^ ,x) ■ 



ax 



bx"^ 
1 



X 



1 + ax + bx"^ ' 1 + ax + bx'^ 



we have 

1 — Ax — /ix^ 



X 



1 + ax + fex^ 1 + ax 
If we now let 



bx"^ 



X 



ax + 6x^ 1 + ax + bx"^ 



Ax — /ix^ ' 



X 



1 



X 



1 + ax' 1 + ax^ 

then [22] we obtain that the Stieltjes matrix for Li is given by 



/A 1 

6 + 10 

6 10 

6 10 

6 1 

6 

V ; ; ; ; ; ; 

We have in fact the following general result [22] : 



\ 



Proposition 11. If L 

necessarily 



{g{x),f{x)) is a Riordan array and P = Sl is tridiagonal, then 



P = Sr 



/ ai 1 

10 

a 1 

6 a 1 

6 a 1 

6 a 



ai 
6i a 
6 
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where 



f{x) = Rev- 



X 



ax 



hx"^ 



and g{x) 



1 

aix — bixf ' 



and vice-versa. 

This leads to the important corollary 

Corollary 12. If L = {g{x), f{x)) is a Riordan array and P = Sl is tridiagonal, with 



Sr 



/ 


ai 


1 











. 






bi 
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1 








. 
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a 


1 





. 












b 


a 
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. 
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a 


1 . 


















b 


a . 




\ 

















(6) 



then L ^ is the coefficient array of the family of orthogonal polynomials Pn{x) where Po{x) 
1, Pi{x) = X — ai, and 



Pn+l{x) = (X - a)pn{x) - bnPn-l{x), 

where bn is the sequence 0,bi,b,b,b, . . .. 



n > 2, 



Proof. By Favard's theorem, it suffices to show that defines a family of polynomials 
{Pn{x)} that obey the above three-term recurrence. Now L is lower-triangular and so 
is the coefficient array of a family of polynomials p„,(x) (with the degree of p„,(x) being n), 
where 



L 



-1 



We have 
Thus 



Sr-L 



-1 



L 



-1 



1 


\ 




( Po{x) 


X 






pi{x) 


x^ 






P2{x) 


x^ 






Ps{x) 




) 






U 


-1 


= L 





We therefore obtain 



/ «! 


1 














bi 


a 


1 














b 


a 


1 














b 


a 


1 














b 


a 


1 














b 


a 



^ = L-^ ■ I ■ {l,x,x\ 
\ 



= L-'-l 



2 s 

Xnn nf*'-' 



/ Po(a;) \ 
Px{x) 
P2ix) 

V '■■ 



( xpo{x) \ 
xpi{x) 

Xp2{x) 
XP2,{x) 
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from which we infer that 

Pi{x) = X — ai, 

and 

Pn+l{x) + apn{x) + hnPn-l{x) = XPn{x) , U > 1, 

or 

Pn+l{x) = {X- a)pn{x) - bnPn-l{x), U > 1. 

□ 

If we now start with a family of orthogonal polynomials {pn{x)}, po{x) = 1, pi{x) = x — ai, 
that for n > 1 obey a three-term recurrence 

Pn+l{x) = (X - a)pn{x) - bnPn-l{x), 

where bn is the sequence 0,bi,b,b,b, . . ., then we can define [23] an associated Riordan array 
L = igix)J{x))hj 

X 1 

fix) = Rev— -— and g{x) = — — . 

1 + ax + ox^ 1 — aix — bixj 

Clearly, L^^ is then the coefficient array of the family of polynomials {pn{x)}- Combining 
these results, we have 

Theorem 13. A Riordan array L = {g{x), f{x)) is the inverse of the coefficient array of a 
family of orthogonal polynomials if and only if its production matrix P = Sl is tri-diagonal. 

Proof. If L has a tri-diagonal production matrix, then by Corollary (12), is the coeffi- 
cient array of orthogonal polynomials. It conversely is the coefficient array of a family 
of orthogonal polynomials, then using the fact they these polynomials obey a three-term 
recurrence and the uniqueness of the Z- and A-sequences, we see using equations (2) and 
(3) along with Proposition 9, that the production matrix is tri-diagonal. □ 

Proposition 14. Let L = {g{x), f{x)) be a Riordan array with tri-diagonal production 
matrix Sl. Then 

mx) = 

where C is the linear functional that defines the associated family of orthogonal polynomials. 
Proof. Let L = {h,j)i,j>o- We have [31] 

n 

X ^ ^ ln,iPi (■^) • 

i=0 

Applying £, we get 

(n \ n n 

^ln,iPi{x) j = ^ln,ijC{Pi{x)) = ^ln,A,0 = L,0 = [x'']g{x). 
i=0 J 1=0 i=0 

□ 
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Thus under the conditions of the proposition, by Theorem 6, g{x) is the g.f. of the 
moment sequence = Hence g{x) has the continued fraction expansion 

1 

9{x) = 



bix 

1 — aix 



1 — ax 



bx 

1 — ax 



1 — ax — ■ ■ 

This can also be estabhshed directly. To see this, we use the 
Lemma 15. Let 

f{x) = Rev — r. 

1 + ax + bx^ 

Then 

f 1 



X 1 — ax — bx'^{f /x) 
Proof. By definition, f{x) is the solution u{x), with n(0) = 0, of 

u 



1 + au + bu^ 

We find 



X. 



1 — ax — \/l — 2ax + (a^ — 46)x^ 

U[X) = 



2bx 

Solving the equation 



v(x) 



1 — Of — bx'^v ' 
we obtain 



^ ^ \ — ax — \/l — 2ax + (a^ — 46)x^ f{x) 



X 



Thus we have 



f{x) 1 X 
-Kev 



X X 1 + ax + bx"^ bx 

I — ax 



bx'^ 

1 — ax 



Now 



1 — aix — bixf 1 — aix — bix'^{f /x) 



□ 
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immediately implies by the above lemma that 



9[x) 



aix 



bix 



1 — ax 



bx^ 



bx 



ax 



1 



We note that the elements of the rows of L ^ can be identified with the coefficients of the 
characteristic polynomials of the successive principal sub-matrices of P. 

Example 16. We consider the Riordan array 

1 X 



ax 



bx"^ 1 + ax + bx"^ 



Then the production matrix (Stieltjes matrix) of the inverse Riordan array ( i^ax+hx'^ ' i+ax+hx'^ ) 
left-multiplied by the /c-th binomial array 



is given by 



P 



and vice-versa. This follows since 
1 X 



ir. 


1 


X ^ 




X 






- kx^ 1 


— kx J 




x' 1 — 


X ) 


+ k 


1 














b 


a + k 


1 














b 


a + k 


1 














b 


a + k 


1 














b 


a + k 


1 














b 


a^k 



ax 



+ bx^ ' 1 + ax + bx^ } \ 1 + /cx ' 1 + /cx 



X 



1 



X 



1 + (a -|- k)x -|- bx^ ' 1 -|- (a -|- k)x + bx^ 



In fact we have the more general result : 

1 + Ax + /ix^ X 



ax + bx^ ' 1 + ax + bx'^ 
1 + Ax + /ix^ 



X 



1 + A;x' 1 

X 



kx 



1 + (a + A;)x + fex^ ' 1 + (a + k)x + 6x^ 
The inverse of this last matrix therefore has production array 
/ a 



+ k-X 


1 














b — fi 


a + k 


1 














b 


a + k 


1 














b 


a + k 


1 














b 


a + k 


1 














b 


a + k 
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4 Exponential Riordan arrays 



The exponential Riordan group [3, 8, 10], is a set of infinite lower-triangular integer matrices, 
where each matrix is defined by a pair of generating functions g{x) = go + gix + g2X^ + ■ ■ ■ 
and f{x) = fix + /2X^ + ■ ■ ■ where go ^ Q and /i ^ 0. In what follows, we shall assume 

90 = fi = 1- 

The associated matrix is the matrix whose i-th column has exponential generating function 
g{x)f{xy/i\ (the first column being indexed by 0). The matrix corresponding to the pair 
/, g is denoted by [g, /]. The group law is given by 

[gJ]-[h,l] = [g{hof)Jof]. 

The identity for this law is / = [1, x] and the inverse of [g, f] is [g, /]~^ = [l/{go f)^ f] where 
/ is the compositional inverse of /. We use the notation eTZ to denote this group. 

If M is the matrix [g, /], and u = (u„)„>o is an integer sequence with exponential generat- 
ing function U (x), then the sequence Mu has exponential generating function g{x)h({f{x)). 
Thus the row sums of the array [g, f] have exponential generating function given by g(x)e^^^^ 
since the sequence 1,1,1,... has exponential generating function e^. 

As an element of the group of exponential Riordan arrays, the Binomial matrix B is 
given by B = [e^,x]. By the above, the exponential generating function of its row sums is 
given by e^e^ = e^^, as expected (e^^ is the e.g.f. of 2"). 

Example 17. We consider the exponential Riordan array [^3^,0;], AQ94587 . This array has 
elements 



/ 1 

















1 


1 














2 


2 


1 











6 


6 


3 


1 








24 


24 


12 


4 


1 





120 


120 


60 


20 


5 


1 



V : : : : : : 

and general term [k < and inverse 

/ 1 

-110 

0-2 1 

0-3 1 

-4 1 

-5 1 

\ ; ■■ ■■ ■■ J 

which is the array [1 — x,x]. In particular, we note that the row sums of the inverse, which 

begin 1, 0, —1, —2, —3, . . . (that is, 1 — n), have e.g.f. (1 — x) exp(x). This sequence is thus 
the binomial transform of the sequence with e.g.f. (1 — x) (which is the sequence starting 
1,-1,0,0,0,...). 
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Example 18. We consider the exponential Riordan array L = [1, jz^]- The general term of 
this matrix may be calculated as follows 



-'- n,k , . \X 



^[^-1(1-^)" 



j=0 



n\fk + n — k — 1 
k\\ n — k 
n\ fn — 1 
k\\n — k 

Thus its row sums, which have e.g.f. exp (yz^); have general term X^fc=o f[ (n-D • This is 
A000262 , the 'number of "sets of lists" : the number of partitions of {1, .., n} into any number 
of lists, where a list means an ordered subset'. Its general term is equal to (n— l)!L„_i(l, —1). 

We will use the following [8, 10], important result concerning matrices that are production 
matrices for exponential Riordan arrays. 

Proposition 19. Let A = {(in,k)nk>o ~ exponential Riordan array and let 

c{y) = Co + Ciy + C2?/^ + . . . , r{y) = Tq + riy + rg?/^ + . . . (7) 
be two formal power series that that 

r{f{x)) = fix) (8) 

= (9) 
9{x) 



Then 



«n+l,0 = ^«!cian,i (10) 

i 

roan,k-i + r[ 5Z ^'(Ci-fc + kri_k+i)an,i (H) 



i>k 



or, assuming = for k < and = for k < 0, 



kl 

i>k-l 



^ ^ i\{ci-k + kri-k+i)an,i. (12) 



Conversely, starting from the sequences defined by (7), the infinite array {cin,k)nk>o defined 
by (12) is an exponential Riordan array. 
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A consequence of this proposition is that the production matrix P = iPi,j)^jyQ for an expo- 
nential Riordan array obtained as in the proposition satisfies [10] 



Pi,j = j{ci-j + jri-j+i) 
Furthermore, the bivariate exponential function 

(j)p{t,z) = y^^pnM^ 

n,k 



0). 



n! 



of the matrix P is given by 
Note in particular that we have 
and 



0p(t,2) = e''{c{z) +tr{z)). 
r{x) = f'{f{x)), 



9{M) ^ 

Example 20. The production matrix of [l, A111596 is given by 



(13) 
(14) 



/o 


1 











. 







-2 


1 








. 







2 


-4 


1 





. 










6 


-6 


1 


. 













12 


-8 


1 . 
















20 


-10 . 




v '■■ 















The row sums of L ^ have e.g.f. exp^j^), and start 1,1,-1,1,1,-19,151,.... This is 
A111884 . This follows since we have g{x) = 1 and so g'{x) = 0, implying that c{x) = 0, and 
fi^) = if^ which gives us f{x) = and f'{x) = jj^- Thus /'(/(x)) = c{x) = (1 - xf. 
Hence the bivariate generating function of P is e^^(l — x)^y, as required. 

Example 21. The exponential Riordan array A = [33^' T3^] ' 



/ 1 

















1 


1 














2 


4 


1 











6 


18 


9 


1 








24 


96 


72 


16 


1 





120 


600 


600 


200 


25 


1 



19 



has general term 



_ru n 
k\\k 



It is closely related to the Laguerre polynomials. Its inverse is the exponential Riordan 
array [j^, with general term (~l)"~'^f[ (^) • This is A021009 , the triangle of coefficients 
of the Laguerre polynomials Ln{x). 

The production matrix of the matrix A^^ = [j^, y+^] given by 



/ 1 1 

13 10 

4 5 1 

9 7 

16 

25 11 

V ; ; ; ; ; ; 







1 

9 







1 



This follows since we have g{x) = and so g'{ 
fix) = if^ and fix) 



X] 



(l-x)- 



r. Then 



c(x) 



9{Kx)) 



(l-x)2 



X 



and /(x) = which yields 



while 



= /'(/») = (l+x)2. 
Thus the bivariate generating function of P is given by 

e'yil + x + il + xfy). 



We note that 



where 



A = exp(S), 



/o 














. 


.. \ 


1 














. 







4 











. 










9 








. 













16 





. 
















25 


. 




v ■■ 
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Example 22. The exponential Riordan array [e^,ln (jz^)], or 



/ 1 


n 


n 


n 


n 


n 


. . ' 


1 


1 











. 




1 


3 


1 








. 




1 


8 


6 


1 





. 




1 


24 


29 


10 


1 


. 




1 


89 


145 


75 


15 


1 . 




V '■■ 















is the coefficient array for the polynomials 

2Fo{-n,x; -1) 

which are an unsigned version of the Charlier polynomials (of order 0) [12, 24, 30]. This is 
A094816 . It is equal to 

' , / 1 

\L — X 

or the product of the binomial array B and the array of (unsigned) Stirling numbers of the 
first kind. The production matrix of the inverse of this matrix is given by 



/ -1 1 
1-21000 

2-3100 

3 -4 1 

4 -5 1 

5 -6 

V I : I I : ! 



which indicates the orthogonal nature of these polynomials. We can prove this as follows. 
We have 

1 



eMn 



X 



Hence g{x) = e ^ and f{x) = 1 — e ^ . We are thus led to the equations 



r(l — e 
c(l-e 



— e 



with solutions r(x) = 1 — x, c(x) = x — 1. Thus the bi-variate generating function for the 
production matrix of the inverse array is 



which is what is required. 



21 



We can infer the following result from the article [23] by Peart and Woan. 

Proposition 23. If L = [g{x), f{x)] is an exponential Riordan array and P = Sl is tridi- 
agonal, then necessarily 



P = Sr 



/ ao 1 

P2 







V '■■ 





1 

^2 









1 

as 







1 

/35 







1 

as 



where {ai}i>Q is an arithmetic sequence with common difference a, is an arithmetic 

sequence with common difference /3, and 



(tto + A/) dx, g{0) 



where f is given by 



and vice-versa. 



f = l + af + Pf\ /(0) = 0, 



We have 



In the above, we note that ao = ai where g{x) is the g.f. of ao = 1, ai, a2, 
the important 

Corollary 24. If L = [g{x), f{x)] is an exponential Riordan array and P = Sl is tridiagonal, 
with 



( 



Sr 





1 














Pi 


ai 


1 














P2 


a2 


1 
















as 


1 


















1 
















as 



v 



\ 



then L ^ is the coefficient array of the family of monic orthogonal polynomials Pn{x) where 
Po{x) = 1, pi{x) = X — ai = X — ao, and 

Pn+i{x) = {x- an)Pn{x) - /3„p„-i(x), n > 0. 

Proof. By Favard's theorem, it suffices to show that defines a family of polynomials 
{pn{x)} that obey the above three-term recurrence. Now L is lower-triangular and so is 
the coefficient array of a family of polynomials Pn{x), where 



/ 1 \ 



X 
2 



X 



V ■■ J 



( Poix) 
pi{x) 

P2{x) 

Ps{x) 

V '■■ 



\ 
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We have 
Thus 

We therefore obtain 



Sr ■ L-^ = L-^ -L-L-^ = L-^ -I-L- L-^ = L"^ ■ /. 



Sl- L ^ ■ 



^ = L-^■J■(l,x,x^...f 



"0 


1 














Pi 




1 


















1 
















"3 


1 


















1 
















as 



v ■■ 

from which we infer 
and 
or 



Po[x) 

P2{x) 



J 



\ 



( 



J 



xpo{x) 
xpi{x) 

Xp2{x) 

xps{x) 



\ 



Pi{x) = x-ao, 

Pn+l{x) + anPn{x) + (3nPn-l{x) = U > I, 

Pn+l{x) = {X- an)Pn{x) - (3nPn-l{x) U > 1. 



□ 



If we now start with a family of orthogonal polynomials {Pn{x)} that obeys a three-term 
recurrence 

Pn+l{x) = {X- an)Pn{x) - (3nPn-l{x) , U > 1, 

with po{x) = 1, pi{x) = X — ao, where {aj}j>o is an arithmetic sequence with common 
difference a and {"T-j^^]^ is an arithmetic sequence with common difference /3, then we can 
define [23] an associated exponential Riordan array L = [g{x), f{x)] by 

/' = 1 + «/ + /?/', /(0) = 0, 



and 



Hg) = j{ao + (3if)dx, (?(0) = 1. 



Clearly, L ^ is then the coefficient array of the family of polynomials {pn{x)}- Gathering 
these results, we have 

Theorem 25. An exponential Riordan array L = [g{x), f{x)] is the inverse of the coefficient 
array of a family of orthogonal polynomials if and only if its production matrix P = Sl is 
tri-diagonal. 

Proposition 26. Let L = [g{x), f{x)] be an exponential Riordan array with tri-diagonal 
production matrix Sl- Then 

nimx) = Cix'') = fxn, 
where C is the linear functional that defines the associated family of orthogonal polynomials. 
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Proof. Let L = {k,j)i,j>o- We have 

n 

Applying C, we get 

(71 \ n n 

4=0 / i=0 i=0 

□ 

Corollary 27. Lei L = [g{x), f{x)] be an exponential Riordan array with tri- diagonal pro- 
duction matrix Sl- Then the moments fin of the associated family of orthogonal polynomials 
are given by the terms of the first column of L. 

Thus under the conditions of the proposition, g{x) is the g.f. of the moment sequence 
£(a;"). Hence g{x) has the continued fraction expansion 

9{x) = ^ 



(3ix 

1 — UqX 



aix 



a2X — 



— a^x — ■ ■ ■ 

When we come to study Hermite polynomials, we shall be working with elements of the 
exponential Appell subgroup. By the exponential Appell subgroup AeTZ of eTZ we understand 
the set of arrays of the form [f{x),x\. 

Let A G AelZ correspond to the sequence {an)n>o, with e.g.f. /(x). Let B G AelZ 
correspond to the sequence (6„), with e.g.f. g{x). Then we have 

1. The row sums of A are the binomial transform of (a„). 

2. The inverse of A is the sequence array for the sequence with e.g.f. j^. 

3. The product AB is the sequence array for the exponential convolution a * b{n) = 
ELo Qo-khn-k with e.g.f. f{x)g{x). 

For instance, the row sums of A = [f{x),x] will have e.g.f. given by 

[fix),x] ■ = f{x)e^ = e^fix) = [e^x] ■ /(x), 
which is the e.g.f. of the binomial transform of (an)- 
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Example 28. We consider the matrix [cosh(x),x], A119467 , with elements 



/ 1 




1 

1 






1 





6 






1 












1 



5 10 1 



v 



The row sums of this matrix have e.g.f. cosh(x) exp(x), which is the e.g.f. of the sequence 
1, 1, 2, 4, 8, 16, . . .. The inverse matrix is [sech(x), x], A119879 , with entries 

\ 



/ 1 




















1 














-1 





1 














-3 





1 








5 





-6 





1 








25 





-10 





1 



V 



The row sums of this matrix have e.g.f. sech(x) exp(x). This is A155585 . 



5 The Hankel transform of an integer sequence 

The Hankel transform of a given sequence A = {oq, Oi, 02, ...} is the sequence of Hankel 
determinants {ho, hi, h2, ■ ■ ■} where hn = |ai+j|"j^Q, i.e 



A = {flnlnGNo h= {hn} 



The Hankel transform of a sequence a„ and its binomial transform are equal. 
In the case that a„ has g.f. g{x) expressible in the form 



(15) 



ao 



1 — aox — 



1 — aix — 



1 — 



^ 

1 — a^x - 



25 



then we have \18] 



n+l-fc 
k 



k=l 



Note that this independent from 



We note that a„, and are in general not integers. Now let H 



determinant of Hankel type with {i,j)-th term fiui+v ■ Let 



Ml 

Vl 



Uk 
Vk 



(16) 



be the 



A„ = H 



1 
1 



n 
n 



K = Hn 



Then we have 



0!r 



A, A, 



/3n 



n-1 



1 ... 
1 ... 

A„_2A,i 
" A2 , ■ 

n— 1 



n — 1 77, 

n — 1 n + 1 



(17) 



6 Legendre polynomials 

We recall that the Legendre polynomials Pn{x) can be defined by 



k=0 



, 2 / -I I \ n~k / -t \ k 

n\ /1+X\ f 1 — X 



Their generating function is given by 



1 



Vl - 2xt + t2 



n=0 



We note that the production matrix of the inverse of the coefficient array of these polynomials 
is given by 

/ 1 ... \ 



I I 

I I 

f I 

I I 

^0 



v 



which corresponds to the fact that the Pn{x) satisfy the following three-term recurrence 

{n + l)P„+i(x) = (2n + l)xPn{x) - nPn-i{x). 
The shifted Legendre polynomials Pn{x) are defined by 

Pn{x)=Pn{2x-l). 
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They satisfy 

Pnix) 



^)"E(:)Cr')<-^)'-E(-r'Ci^-)( 

fc=0 ^ / ^ / k=0 ^ ^ ^ 



'2k 
k 



X 



Their coefficient array begins 



1 

















-1 


2 














1 


-6 


6 











-1 


12 


-30 


20 








1 


-20 


90 


-140 


70 





-1 


30 


-210 


560 


-630 


252 



V ■■ 

and so the ffist few terms begin 



1, 2x - 1, 6x^ - 6x + 1, 20x^ - 30x^ + 12x - 1, . . . 



We clearly have 



n=0 



^l-2{2x-l)t + t^ 

7 Legendre polynomials as moments 

Our goal in this section is to represent the Legendre polynomials as the ffist column of a 
Riordan array whose production matrix is tri-diagonal. We ffist of all consider the so-called 
shifted Legendre polynomials. We have 

Proposition 29. The inverse L of the Riordan array 

1 + r(l — r)x^ X 



1 + (2r — l)x + r(r — l)x'^ 1 + (2r — l)x + r(r — l)x'^ 

has as its first column the shifted Legendre polynomials Pnij). The production matrix of L 
is tri-diagonal. 

Proof. Indeed, standard Riordan array techniques show that we have 

.2 - \ -1 



1 + r( 1 — r)x 



X 



1 + (2r — l)x + r(r — l)x'^ ' 1 + (2r — l)x + r(r — l)x'^ 

1 1 - (2r - l)x - v/F^2(27^^T)^^M^ 



'1 - 2{2r - l)x + x^ 2r{r-l)x 
This establishes the ffist part. Now using equations (1), with 



fix) 



1 - (2r - l)x - v^l - 2(2r - l)x + 



x^ 



2r(r — l)x 



X 



1 + (2r — l)x + r{r — l)x^ 



27 



and 



1 



1 - 2(2r - l)x + x^ 
we find that 

Z{x) = (2r - 1) + 2rx{r - 1), A{x) = 1 + x(2r - 1) + xV(r - 1). 
Hence by Proposition 9 the production matrix P = 5^ of L is given by 



/ 2r - 1 1 

2r(r-l) 2r-l 1 

r(r-l) 2r-l 1 

r(r - 1) 2r - 1 





V : : : 








1 
r(r — 1) 2r — 1 1 

r(r-l) 2r-l 



□ 

Corollary 30. The shifted Legendre polynomials are moments of the family of orthogonal 
polynomials whose coefficient array is given by 



1 + r(l — r)x 



X 



1 + (2r — l)x + r(r — ' 1 + (2r — l)x + r(r — 
Proof. This follows from the above result and Proposition 14. 

Proposition 31. The Hankel transform of the sequence Pn{r) is given by 2"(r(r — 1 
Proof. From the above, the g.f. of -Pn(r) is given by 

1 



□ 



1 - (2r - l)x 



2r(r — l)x 



1 - (2r - l)x 



r{r — \)x 



1 - (2r - l)x 



r{r — l)x 



The result now follows from Equation (16). 
We note that 



□ 



Pn{r) 



X 



2^r(r-l)+2r-l 
-2^r(r-l)+2r-l ^/-X^ + 2(2r - l)x - 1 



dx 



gives an explicit moment representation for Pn{r). 

Turning now to the Legendre polynomials Pn{x), we have the following result. 
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Proposition 32. The inverse L of the Riordan array 

1 I l-r^^2 



-X 



X 



1 + rx + ^^^2^x2 ' 1 + rx + ^^^j^x^ 

/ifls as its first column the Legendre polynomials Pn{r). The production matrix of h is tri- 
diagonal. 

Proof. We have 

' 1 + '-^x^ X 



1 + rx + ^^r^x^ ' 1 + rx + ^^^TT-x'^ 



2(1 — rx — Vl — 2rx + x^ 



Vl - 2rx + x2' x(r2 - 1) 

This proves the first assertion. Using equations (9) again, we obtain the following tri-diagonal 
matrix as the production matrix of L: 



/ r 1 



2 





H-1 



0^ 
c ^ 






10 



1 
r 1 



r 



□ 



Corollary 33. 



1 + ^x^ 



X 



1 + rx + ^^^Y^x'^ ' 1 + rx + ^^^j^x"^ 

is the coefficient array of a set of orthogonal polynomials for which the Legendre polynomials 
are moments. 

Proposition 34. The Hankel transform of Pn{r) is given by 

(r^-l)m 

Proof. From the above, we obtain that the g.f. of P„(r) can be expressed as 



1 — rx — 



1 — rx 



1 — rx 



The result now follows from Equation (16). 



□ 
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We end this section by noting that 



PJr) 



TT 



X 



-\/r 



-1 \/—x^ + 2rx — 1 



dx 



gives an exphcit moment representation for P„(r) 



8 Hermite polynomials 

The Hermite polynomials may be defined as 



I - 1 



Hjx) 



E 

k=0 



-l)'^(2x) 



k\{n-2k)\ 



The generating function for Hn{x) is given by 



^2xt-t^ 



f 

„ n! 



n=0 



The unitary Hermite polynomials (also called normalized Hermite polynomials) are given by 



He^{x) = 2-^H^{V2x) = ^ J"' 



k=0 



Their generating function is given by 



f2 

Tt— — 

e 2 



n=0 



We note that the coefficient array of ife„ is a proper exponential Riordan array, equal to 



e 2 X 



This array A066325 begins 



/ 1 



-1 

3 





1 


-3 


15 





1 




-6 






1 












1 



-10 1 



30 



It is the aeration of the alternating sign version of the Bessel coefficient array A001497 . The 
inverse of this matrix has production matrix 



/ 1 

10 10 

2 1 

3 1 

4 1 

5 



V : : : : : : •• / 

which corresponds to the fact that we have the following three-term recurrence for ife„: 



Hen+i{x) = xHen{x) - nHen-i{x). 



9 Hermite polynomials as moments 



Proposition 35. The proper exponential Riordan array 



has as first column the unitary Hermite polynomials Hcnir). This array has a tri-diagonal 
production array. 

Proof. The first column of L has generating function e''^ 2 , from which the first assertion 
follows. We now use equations (13) and (14) to calculate the production matrix P = Sl- 

We have /(x) = x, so that f{x) = x and f'{f{x)) = 1 = r(x). Also, g{x) = e^^ 2 implies 
that g'{x) = g{x){r — x), and hence 



c{x) 



9'{Rx)) 



r — X. 



Thus the production array of L is indeed tri-diagonal, beginning 



r 1 

— 1 r 

-2 
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r 1 

-3 r 1 



-4 r 1 
-5 r 



v ■■ 



\ 



□ 
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We note that L starts 

/ 1 

r 

r2-l 
r(r2 - 3) 
- 6r2 + 3 



r(r^ - lOr^ + 15) 5(r4 





1 

2r 1 

3(r2 - 1) 3r 1 

4r(r2 - 3) Q{r^ - 1) 4r 1 

6r2 + 3) 10r(r2 - 3) 10{r^ - 1) 5r 








1 



Thus 



is the coefficient array of a set of orthogonal polynomials which have as moments the unitary 
Hermite polynomials. These new orthogonal polynomials satisfy the three-term recurrence 

with 9)0 = 1, S^i = X — r. 

Proposition 36. The Hankel transform of the sequence Hcnij) is given by (—1)^ ^ ) Y\!k=o ^' 
Proof. By the above, the g.f. of Hn{r) is given by 

1 



1 — rx 



2x' 



1 — rx 



3x^ 



1 — rx + 



The result now follows from Equation (16). 

Turning now to the Hermite polynomials Hn{x), we have the following result. 
Proposition 37. The proper exponential Riordan array 



□ 



Jlrx—x^ 



X 



has as first column the Hermite polynomials Hn{r). This array has a tri-diagonal production 
array. 

Proof. The ffist column of L has generating function e^™~^^ , from which the ffist assertion 
follows. Using equations (13) and (14) we find that the production array of L is indeed 
tri-diagonal, beginning 



/ 2r 


1 














-2 


2r 


1 














-4 


2r 


1 














-6 


2r 


1 














-8 


2r 


1 














-10 


2r 
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□ 



We note that L starts 
1 

2r 

2(2r2 - 1) 

4r(2r2 - 3) 
4(4^3 _ i2r^ ^ 3) 




1 

4r 

6(2r2 - 1) 
16r(2r2 - 3) 





1 

6r 

12(2r2 - 1) 







1 











1 



8r(4r4 - 20r^ + 15) 20(4r4 - Ur"^ + 3) 40r{2r^ - 3) 20(2r2 - 1) lOr 1 



V 

Thus 



-1 



X 



is the coefficient array of a set of orthogonal polynomials which have as moments the Hermite 
polynomials. These new orthogonal polynomials satisfy the three-term recurrence 

Sjn+i{x) = (x - 2r)f)„(x) + 2nS)n^i{x), 

with Sjq = 1, Sji = X — 2r. 

Proposition 38. The Hankel transform of the sequence Hn{r) is given by (—1)^ ^ ) YYk=o 2'^^' 
Proof. By the above, the g.f. of Hn{r) is given by 

1 



1 — 2rx 



2x' 



1 — 2rx + 



4x^ 



1 — 2rx 



6x^ 



The result now follows from Equation (16). 



□ 



10 Acknowledgements 

The author would like to thank an anonymous referee for their careful reading and cogent 
suggestions which have hopefully led to a clearer paper. 



11 Appendix - The Stieltjes transform of a measure 

The Stieltjes transform of a measure on M is a function defined on C \ M by 

1 



z-t 



fi{t). 
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If / is a bounded continuous function on M, we have 

/ = — lim / f{x)QGfj_{x + iy)dx. 

If fi has compact support, then is holomorphic at infinity and for large z, 

oo 



n=0 

where a„ = J^t"'n{t) are the moments of the measure. If fi{t) = dip(t) = ip'{t)dt then 
(Stieltjes- Perron) 

1 /•* 

ip{t) — ip{tQ) = hm / '^GJx + iy)dx. 

71" y^o+ J to 

If now g{x) is the generating function of a sequence a^, with g{x) = J2'^=o '^n^^", then we can 
define 



^ ^ n=0 



By this means, under the right circumstances we can retrieve the density function for the 



measure that defines the elements a„ as moments. 



Example 39. We let g{z) = ^ be the g.f. of the Catalan numbers. Then 



Then 



V2J y/x^ + y^^x^ - 8x + 7/2 + 16 - x2 + 4x - 
+ iy) = 



4:^/x^~+y^ 



and so we obtain 



X j \/x2 + ?/2 a/x^ — 8x + + 16 — + 4x — ?/2 
'(/''(x) = lim < ^ 



_ 1 y/x{A - x) 
27r X 
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